Footprints of sticky motion in the phase space of 
higher dimensional nonintegrable symplectic 
systems 



C. Manchein 1 , M. W. Beims 1 2 and J. M. Rost 2 

1 Departamento de Ffsica, Universidade Federal do Parana, 81531-980 Curitiba, PR, 
. Brazil 
£NJ | 2 Max Planck Institute for the Physics of Complex Systems, Nothnitzer Strasse 38, 

£>^' D-01187 Dresden, Germany 

E-mail: mbeims@fisica.ufpr.br 

. Abstract. "Sticky" motion in mixed phase space of conservative systems is difficult 

to detect and to characterize, in particular for high dimensional phase spaces. Its effect 
f""^ ' on quasi-regular motion is quantified here with four different measures, related to the 

. distribution of the finite time Lyapunov exponents. We study systematically standard 

maps from the uncoupled two-dimensional case up to coupled maps of dimension 20. 
We find that sticky motion in all unstable directions above a threshold Kd of the 
nonlinearity parameter K for the high dimensional cases d = 10, 20. Moreover, as 
K increases we can clearly identify the transition from quasircgular to totally chaotic 
(f) . motion which occurs simultaneously in all unstable directions. The results show that all 

^ [ four statistical measures sensitively probe sticky motion in high dimensional systems. 
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1. Introduction 



The dynamics of conservative systems consists of different types of motion, from regular 
over quasi-regular to chaotic, depending on the perturbation parameter. Usually, in 
d = 2 dimensional conservative systems regular islands (Kolmogorov-Arnold-Moser 
KAM tori) can break as the perturbation parameter increases, and chaotic trajectories 
may coexist with regular islands [jQ. This regime is called quasiregular. Chaotic 
trajectories may be trapped around stable islands creating "sticky" motion (for a review 
see [2]). In such cases the dynamics in phase space can be very structured and the motion 
is non-ergodic. Analytical results for physical quantities in the quasiregular regime 
are rare and it is very desirable to obtain tools which quantify degrees of ergodicity 
and hyperbolicity, as for example the angle between the stable and unstable manifolds 
obtained via the covariant Lyapunov vectors [3l H] . Even for apparently totally chaotic 
systems, tiny regular islands may appear which induce sticky motion as exemplified 
in Fig. with the phase space of the d = 2 dimensional kicked rotor [Eqs. ([T] for 
the nonlinearity parameter K = 5.79. The dynamics looks totally chaotic - yet, 
tiny regular islands appear near the points (x,p) ~ (4.6, 2.9), (1.7, —2.9) (see boxes), 
if enough initial conditions are sampled. Such islands become visible with sufficient 
magnification (Figs. [Tb-c). Hence, the question arises how these tiny regular islands can 
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Figure 1. Phase space for K = 5.79. (b) and (c) are magnifications from the boxes 
shown in (a). 



be identified with not more information than provided by Fig. [TJi. In higher dimensions, 
the information equivalent to Fig. [Th, must be cast in a different form since Poincare 
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Surfaces of Section (PSS) cannot be constructed. In fact, for d > 2 dimensional 
conservative systems, due to Arnold diffusion, the KAM surfaces do not divide the 
whole phase space into a set of closed volumes. Although the effect of stickiness of a N- 
dimensional invariant tori has been analyzed for the asymptotic decay of correlations and 
anomalous diffusion [5], to the best of our knowledge, the behavior of sticky trajectories 
in higher-dimensional systems has never been quantified and explored systematically. 

In the following we present extensive numerical simulations to detect and quantify 
sticky motion in the phase space of quasiregular and chaotic conservative systems with 
up to 20 dimensions. To fulfill this task we analyze small changes in the finite time 
Lyapunov exponents (FTLE) distribution using four quantities: the variance a, the 
skewness k$, the kurtosis k 4 and Pa, the normalized number of occurencies of the 
most probable finite time Lyapunov exponent. The first three quantities were analyzed 
in the two-dimensional standard map for the distribution of the finite time stability 
exponents to detect islands [5] and to study the fluctuations of the mean FTLE for 
specific trajectories and the corresponding Kolmogorov like entropy [TJ. 

The FTLEs are the average of local divergence of trajectories and for ergodic 
systems they converge to the LE obtained in the limit t — > oo, which is independent of 
initial conditions [8]. In such cases the FTLE distribution has Gaussian form (normal 
distribution) [HI HOI EI]- However, in quasiregular systems the chaotic trajectories may 
visit the neighborhood of regular islands, being trapped there for a while, and the FTLEs 
do not converge quite well. Such small changes in the FTLEs distribution can be used for 
many purposes: to obtain information about algebraic or exponential stretching in phase 
space [T2J [T3] , to quantify ergodicity [HI EE] , to study stretching processes in chaotic 
flows [16] , to detect the bimodal [17] or multimodal [18] behavior of the local FTLEs, to 
analyze the time correlations [191 I2Q1 121] and Poincare recurrence [21] in coupled maps, 
to detect the sticky motion [221 ES] and the arise of regular islands [23] for interacting 
particles in one- and two-dimensional billiards via Pa and, as we will see here, to detect 
and quantify sticky motion in higher- dimensional nonintegrable conservative systems. 
In addition, FTLEs have been used in dissipative systems to distinguish between strange 
chaotic and nonchaotic attractors (see [251 126], (27J [28] for some examples). In these cases 
the trajectory switches between long times on the original periodic attractor (negative 
FTLE) and smaller times on the chaotic saddle with positive FTLEs. This is quite the 
opposite of what is observed in conservative systems, where the trajectory stays longer 
times along the chaotic trajectory and smaller times near the regular island, which than 
reduces the FTLE. 

Different from the large literature related to the analysis of sticky motion via the 
distribution of the local FTLEs for some specific trajectories, in this work we analyze 
the FTLEs distribution after n = 10 7 iteration of the map and over 10 4 initial conditions 
distributed uniformly in phase space. Our investigations show that even for very large 
times the FTLEs still depend on initial conditions when sticky motion is present. The 
paper is organized as follows. In Section [2] we describe the properties of the variance a 
and the distribution of the most probable Lyapunov exponent Pa for a two-dimensional 
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conservative system (the standard map) and compare them with higher cummulants (^3 
and k 4 ). In Section [3] we couple N standard maps and by increasing N = 2,3, 5, 10 we 
analyze carefully the changes observed in the the quantities Pa, k% and k^. Section |4] 
summarizes the main results of the paper. 

2. The standard map and the FTLE distributions 

The standard map is an ideal model for a conservative chaotic system to study the 
behavior of FTLE distributions, especially if one is interested in an extension to higher 
dimensions. The original d = 2 dimensional map M was proposed by Chirikov [29J and 
is defined by 

J Pn+i = p n + K sin(x n ) mod 2vr, 

' I x n+1 = x n + p n+ i mod 2ix, 

where p n and x n are the dynamical variables of the system after the nth iteration and 
K is the nonlinear parameter. This map locally approximates more general nonlinear 
mappings pQ. The FTLE spectrum is determined for a symplectic system using the 
standard procedure [301 EI], consequently the sum of all FTLEs is zero. This implies 
that the FTLEs come in pairs ±A^ and that the number of zero exponents is even [i 
is for the zth unstable direction). Hence, it is sufficient to analyze the An^ > 0. The 
FTLEs can be defined by 

A^(x) = ^ln||M(x,n)e 4 (x)||, (2) 

with the discrete iterations n. M(x,n) = J(x n _i) • • • J(xi) J(x ), and J(x n ) (n = 
0, . . . , n — 1) are the Jacobian matrices determined at each iteration n. ej(x) represents 
the unit vector in the ith unstable direction. These vectors are re-orthonormalized at 
each iteration and the FTLE spectrum is obtained using Eq. (J2J). The usual LE is 
obtained from A 1 = lim n ^oo A^(x) which is independent of the initial conditions. 

Figure^ shows the FTLEs distributions P(A n , K) as a function of the nonlinearity 
parameter K for 400 initial conditions and up to n = 10 7 iterations. For each initial 
condition we reckon one FTLE A n , which is the time average of all local FTLE along the 
trajectory. In Fig. [2^ the color changes from light to dark (white over yellow and blue to 
black) with increasing P(A n ,K). The gray points below the main curve are related to 
FTLEs from chaotic trajectories which were trapped for a while close to regular islands. 
For these trapped chaotic trajectories (or "sticky" motion) the local FTLE decreases 
and consequently the final A n also decreases. In such cases the FTLE distribution is not 
Gaussian anymore but contains interesting information about the dynamics. Figure [2^ 
also reveals that for some specific values of K the amount of sticky trajectories changes, 
which will be quantified and characterized next. Note, that such information is not 
provided by the mean FTLE shown in Fig. [2b. 
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Figure 2. (a) (Color online) Finite-time distribution of the positive Lyapunov 
exponent P(A„,K) calculated over 400 trajectories up to n = 10 7 iterations. With 
increasing P(A n ,K) the color changes from light to dark (white over yellow and blue 
to black), (b) The corresponding mean value (A„). 



= 0, (3) 

a„=aS 



2.1. Quantitative characterization of the distributions 

Clearly, anomalies in the FTLE distribution reflect sticky motion - however, the entire 
distribution is not a convenient tool to identify and measure sticky motion. Instead, we 
will analyze the following four measures: 

a. ) The number of occurencies of the most probable FTLEs -P(A^, K) = P\(K), defined 

by [22J 1231 [23] 

dP(A n ,K) 
dk n 

is the value of the maximum of P(A n ,K) for each K (see Fig. [2]). P^{K) = 1 
means that for all initial conditions the FTLEs are equal (within the precision 
10~ 3 ). Therefore we expect P\(K) = 1 for a totally chaotic system and large times. 
For cases of sticky motion appears or if the normal distribution did not converge 
to the given precision, Pa(K) < 1. 

b. ) The variance defined by 

a = <(A„-(A n ,» 2 ), (4) 

which is the second cummulant of the FTLE distribution. It should increase for 
sticky motion. We will analyze the properties of the relative variance o = a / (A„) 2 
which is more appropriate for higher-dimensional systems. In such cases it is 
possible to detect small differences relative to each unstable direction. 
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c. ) The skewness is defined by 

((An ~ (An)) 3 ) „ 

K 3 = • (5) 

It is the third cummulant of the FTLE distribution and detects the asymmetry of 
the distribution around its mean value. For K3 = we have the regular distribution. 
Since sticky motion usually reduces the FTLE, the asymmetry of the distribution 
leads to k 3 < 0. 

d. ) The kurtosis is defined by 

((An - (A»)) 4 ) 

«4 = ~n 3, (6) 

and detects the shape of the distribution, where K4 > indicates that the 
distribution is flatter than the regular distribution (k 4 = 0) while < reveals a 
sharper distribution than the normal one. 

These four quantities will be used to characterize the degree of sticky motion in the 
transition from low- to higher-dimensional systems. Altogether they should give all 
relevant informations which can be extracted from the distributions. 



2.2. The two-dimensional case (d = 2) 

To compare with the higher dimensional cases to be discussed later we summarize here 
d = 2 results for the FTLEs distributions of the standard map in a slightly different K 
regime than considered by [6]. For all simulations we have used 201 values of K in the 
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Figure 3. (Color online) Deviations from the mean Lyapunov exponent in the interval 
K = (7.0,9.0) for the standard map measured with (a) P\(K), (b) <r(K), (c) n^{K) 
and (d) k 4 (K), see Eq. (pTTj) . 

interval K = (7.0, 9.0) (stepsize AK = 0.01). As can be seen from Fig. [3^ Pa has many 
minima in the interval K = (7.0, 9.0). For all minima of Pa a certain amount of sticky 
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motion is expected. This is indeed true, as we have checked constructing the PSS (not 
shown) for some values of K where minima in Pa appear and for all of them we found 
sticky motion. Comparing Pa with Fig. [3b we see that a has maxima only for values 
of K in the interval (7.15,7.55) and close to 8.12, i.e., the variance does not detect 
sticky motion but gives important additional information, as we will see later. However, 
skewness (Fig. [3b ) and kurtosis (Fig. |3H) detect sticky motion as identified by Pa in 
(Fig. 12k). Both parameters behave as expected. In Fig. [3b we plotted the skewness 
«3 as a function of K. For almost all minima of Pa we observe that the skewness 
(Fig. [3b) is negative. This implies that the sticky motion induces a left asymmetry in 
the FTLEs distribution. In other words, the tail on the left side of the distribution 
is larger than the tail on the right side. For almost all minima of Pa the kurtosis k 4 
(Fig. |3H) has a maximum. This is expected since a small Pa implies widely distributed 
A n and consequently a relatively flat distribution (large K4). 

Obviously, in the two-dimensional case Pa, and the cummulants K3 and k 4 are 
sensitive and consistent detectors of sticky phase space trajectories, with the latter two 
being more sensitive, while Pa is less sensitive, but gives additional information (see, e.g., 
the interval K = (7.2, 7.5)), which will become more obvious in the higher dimensional 
cases. 



3. Coupled standard maps 

Coupled standard maps are convenient systems to investigate sticky motion in higher 
dimensional phase space systematically. The equations for N coupled standard maps 
read 

q J Pn+l = Pn ] + K Sin(xl fc+1) - X { n ] ), 

b ■ \ -(*) _ T (*)+ n (*) [n 

where k = 1, . . . , N. We consider periodic boundary conditions p( N+1 > = p^ l \ x^ N+1 ' = 
x^ and unidirectional next nearest neighbor coupling. In Eq. ([7]) we must use modulus 
2ir. The model is invariant under a translation of the x-coordinate. For N = 2 it 
is Hamiltonian, symplectic and conserves the total momentum. For N > 2 it is just 
symplectic. The Lyapunov spectrum always possesses two null exponents for N > 1. 



3.1. N = 2 (d = 4): The separable case 

In this case the system Eq. (J7J) has one positive, one negative and two zero FTLEs. 
We show results related to the positive FTLE in Fig. HI Clearly, the mean value (A n ) 
increases with K (Fig. 2^). For N = 2 the equations of 5* (0) decouple in center of mass 
and relative coordinates. The map in the relative coordinate is the standard map but 
with parameter 2K. Therefore, the curve in Fig. H^i should match those obtained for 
the standard map with parameter values IK = (14, 18). To check this we plotted (A n ) 
in Fig. (dashed line) for the standard map with values 2K = (14, 18). Both curves 
are identical. For all minima of Pa (Fig. |4b) the skewness K3 (Fig. HU) and the kurtosis 
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Figure 4. (Color online) Comparison of (a) (A„) (the curve N = 1 shows results for 
the standard map calculated at 2K = (14,18), see text) with (b) P\ n (K), (c) cr(K), 
(d) k s (K) and (c) k 4 (K) in the interval K = (7.0, 9.0) and for N = 2. 



Ki (Fig. Hk) have also minima and maxima, respectively. At the corresponding values 
of K many trapped trajectories are expected. Since this case is the standard map with 
larger nonlinearity parameter (2K), we do not discuss further details. 

3.2. N = 3 {d = 6): The largest FTLE is most sensitive to sticky motion 

The case iV = 3 has six FTLEs, two positive, two negative and two zero. The analysis 
will be done for the distributions of the two positive FTLEs which are labeled such 
that An'' > An 2 ' 1 . For N = 3 the interval (7.0,9.0) does not exhibit structures in 
the FTLE distributions pointing to sticky trajectories. Hence, we analyze the interval 
K = (0.01, 2.0) which is close to the non-chaotic case K = 0. Figure shows the mean 
values of the first and second FTLEs as a function of K. Both FTLEs increase with K, 
but An' 1 increases faster and has a minimum close to K ~ 0.2. Close to this value A^ 
starts to be finite. Figures E] and [7] show the quantities (a) P^(K), (b) (d) k^\k) 
and (d) 4\K) as a function of K for An ^ and An 2 \ respectively. By comparing 
these plots we see that both distributions detect the sticky motion for 0.1 < K < 0.7. 
In this region the minimum value of k$(K) is ~ —50 for i — 1 and —30 for % = 2. The 
maximum value of k^(K) is 4000 for i = 1 and ~ 2000 for i — 2. For fT < 0.1 we 
observe in Fig. [7| that P^\K) = 1, meaning that all initial conditions leave to the same 

(2) (2) (2) (2) 

0.0. However, k 3 (K) and k\ (K) show that the distribution around A n ~ 0.0 
is not Gaussian. By changing K = 0.01 — > 0.1 the distribution changes slightly from 
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Figure 5. (Color online) Mean values for the two positive FTLEs A« > ^-n' for 
N = 3. 
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Figure 6. (Color online) Comparison of (a) pP(K) with (b) a^{K), (c) k£\k) 
and (d) k£\K) in the interval if = (0.01, 2.0) and JV = 3. 

f2l (2^ 

the right asymmetric form [k 3 (0.01) ~ 1] to the left asymmetric form [k 3 (0.1) ~ —3]. 

(2) 

The kurtosis is always flat (k\ > —3) but approaches a regular distribution [k = 0) 
when K ~ 0.05. Similar behavior is observed in Fig. [(2 for k^\k) and k± (K). For the 
relative variance is shown to be relevant only at very low values of K (for both FTLEs). 

For K > 0.7 the distribution related to An detects the sticky motion at K ~ 
0.83,1.40 and K ~ 1.32, while for A^ 2) it detects K ~ 1.40 and K ~ 0.83 (but 
only k^\k)). This shows the intuitively expected result that when the nonlinearity 
parameter increases the most unstable dimension (related to A^) is more sensitive to 
the sticky motion. 

3.3. N = 5(d = 10): The smallest FTLE defines a critical nonlinearity Kd 

In Fig. [8] the mean values for the positive FTLEs A n 1} > A^ > A { n ] > An ] are shown 
as a function of K. For smaller values of K < 0.15, < A^ > and < A^ > are almost 
zero, while < A„ > increases initially very fast for K = 0.01 — > 0.05 and then slowly 
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Figure 7. (Color online) Comparison of (a) P}?\K) with (b) a (2 \K), (c) K^\K) 
and (d) k^{K) in the interval if = (0.01, 2.0) and JV = 3. 




Figure 8. (Color online) Mean values for the four positive FTLEs > > 
Al 3) > A< 4) for N = 5 . 



for K = 0.05 — > 0.15. For K above 0.2 all mean FTLEs increase linearly with different 
slopes size-wise ordered as the FTLEs themselves. 

Figure H2 shows that for K > 0.3 the distributions from all FTLEs approach 
the normal (Gaussian) form expected for a totally chaotic regime in contrast to the 
cases N = 1,2,3 discussed before, where sticky motion also occurs for higher values of 
K. Therefore, there is a clear transition from the quasiregular to the totally chaotic 
dynamics at K c ~ 0.3 simultaneously in all unstable directions. On the other hand, all 
quantities detect the sticky motion for K < 0.25. Interestingly, above a critical value 
Kd ~ 0.13, the four indicators of sticky motion have a similar qualitative behavior for all 
FTLEs, i. e., independently of the unstable directions. For example, the distributions 
for all FTLEs have a small left asymmetry («3 ~ —0.75) with no flatness (k^ ~ 0) at 
K = 0.13, and a left asymmetry (ac 3 ~ —1.50) with flatness (k 4 ~ 3.0) at K = 0.2. The 
critical value ~ 0.13 is the point where P^\k) detects the sticky motion from the 
distribution of the last (the smallest) FTLE [see Fig. |9K|- The value of Kd can also be 
determined approximately from the onset of linear .fT- dependence of the mean FTLEs 
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Figure 9. (Color online) (a) P^(K), (b) a^{K), (c) K^{K) and (d) (K) in the 
interval K = (0.01, 1.0), N = 5 and for all i = 1,2,3,4. 



in Fig. El 

These results show that the sticky motion induces differences between the 
distributions of only for K < K d . Since sticky motion is a consequence of 
regular structures in phase space, this suggests that there could exist disconnected 
invariant regions in phase-space. In the interval < K < K c on the other hand, the 
sticky motion affects all unstable dimensions by the same amount leading to similar 
distributions for different FTLEs. Finally, in the totally chaotic regime K > K c all 
unstable directions are affected proportionally to each other. To explain this better we 
show in Fig. flOl the time evolution of h$ (i = 1,2,3,4) for one exemplary trajectory 
in the three distinct regions of K. Figure ITUT a) shows the behavior of the FTLEs 
for K = 0.6 > K c and no sticky motion is observed since no FTLE goes to zero 
for some shorter time intervals. For the intermediate case < K = 0.15 < K c 
we see in Fig. ITUT b) that for times around 0.5 x 10 6 (see arrow) all four FTLEs are 
zero, so the sticky motion occurs in all unstable directions and a common behavior is 
observed. Physically what happens is that the linear momentum of each site becomes 
approximately constant during sticky motion. This is clearly shown in Fig JTTT a). where 
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Pn (j = 1)2,3,4,5) is plotted as a function of n. At the sticky time n ~ 0.5 x 10 6 
(see black box), pn is constant for all sites. This agrees with the conjecture made 




Figure 10. (Color online) Time evolution of A„ (i = 1,2,3,4) for one exemplary 
trajectory with N = 5 and (a) K = 0.6 > K c , (b) K d < K = 0.15 < K c , (c) 
K = 0.05 < K d . 



by Froeschle [32] and observed later in coupled maps by [20], i. e., that a conservative 
system with N > 3 has iV or one (the energy) constants of motion. In other words, 
following this conjecture when the sticky motion occurs all FTLEs should temporarily be 
very close to zero, as already observed. However, this changes when the region K < 
is analyzed, as can be observed in Fig. [TOTc) for times around 0.5 x 10 6 and 0.85 x 10 6 
(see arrows). Only is larger and the remaining three FTLEs are very close to zero. 
In this case, at the sticky time n ~ 0.5 x 10 6 (see black box), the momenta pn\p^ and 
Pn^ are not conserved anymore. See Fig. [TTTb). 



3.4. N = 10(d = 20) 

The case N = 10 has nine positive FTLEs. In this section we analyze again the interval 
K = (0.01, 1.0) but with 10 3 initial conditions only. In Fig. [12] the mean values for the 
FTLEs are plotted as a function of K. For K > 0.2 all mean FTLEs increase linearly 
and the largest ones increase faster. For K < 0.2 the K dependence is different for 
different FTLEs. For K < 0.1 almost all (but two 2 = 1,2) are very close to zero. 

Figure d3] compares (a) Pjp(K), (b) a®(K), (c) 4\ K ) and ( d ) ^t\ K ) m the 
interval K = (0.01, 1.0) for all positive FTLEs (z = 1 9). For K > K c ~ 0.28 all 
quantities indicate that the distributions from all FTLEs approach the normal form 
expected for a totally chaotic regime. Therefore, again we have a clear transition 
from quasi-regular to totally chaotic dynamics at K c which occurs simultaneously in 
all unstable directions. For K < K c all quantities detect the sticky motion. Also here 
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Figure 11. (Color online) Time evolution of Pn (j — 1,2,3,4,5) for the trajectory 
from fig EH (a) K d < K = 0.15 < K c , (b) K = 0.05 < K d . Black boxes show the sticky 
region. 
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Figure 12. (Color online) Mean values for the nine positive FTLEs A^ 1 ' > A%> > 
... > AL 9) , for N = 10. 



,(2) 



we observe that for K above the critical value ~ 0.09 the four quantities behave 
similarly, independent of the FTLEs. The distributions for all FTLEs have a small left 
asymmetry (^3 ~ —0.75) with no flatness (k^ ~ 0) at K = 0.13, and a left asymmetry 
(k 3 ~ —1.50) with flatness (k 4 ~ 3.0) at K = 0.2. The critical value K d ~ 0.09 is 
again exactly the point where P^\k) detects the sticky motion from the distribution 
of the last (the smallest) FTLE [see Fig. [9]^a)]. It can be observed that the distribution 
of the largest FTLE is the last (when K increases from 0.01) to "join" the collective 
behavior for Kd > 0.09 [please see Fig. ITBTbWd)]. The value Kd can also be determined 
approximately from the region in Fig. [T5] where the K dependence of the mean FTLEs 
is not linear. 

Figure HH shows the time evolution of A„ (i = 1 — > 9) for one exemplary trajectory 
in the three distinct regions of K. In Fig. [T47 a) with K = 0.6 > K c no sticky 
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Figure 13. (Color online) Quantities (a) P ( ^\k), (b) cr {i) {K), (c) t$(K) and (d) 
k^{K) as a function of K = (0.01, 1.0), for i = 1 ->■ 9 and TV = 10. 



motion is observed since no FTLE goes to zero for some shorter time intervals. For 
the intermediate case < K = 0.15 < K c we see in Fig. [T4T b) the two different 
scenarios as conjectured. For times around 0.22 x 10 6 (see arrow) all nine FTLEs are 
close to zero indicating that the sticky motion occurs in all unstable directions. For 
times around 0.8 x 10 6 (see arrow) on the other hand, eight FTLEs go to zero but the 
largest one remains positive. This behavior is observed even more clearly in the region 
K < K d , see Fig. 23(c). For times in the interval 0.12 x 10 6 to 0.18 x 10 6 only is 
large and the remaining eight FTLEs are very close to zero. 

4. Conclusions 

It is well known that the dynamics of conservative systems consists of qualitatively 
different types of motion, from regular over quasi-regular to chaotic, depending on 
the perturbation parameter. It is difficult to identify these types of motion in high- 
dimensional systems, especially in the quasi-regular regime where the sticky motion can 
suppress Arnold diffusion and distinct stochastic regions with different properties may 
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Figure 14. (Color online) Time evolution of A„ (i = 1 — > 9) for one exemplary 
trajectory with N = 10 and (a) A" = 0.6 > K c , (b) A' d < K = 0.15 < A c , (c) 
A' = 0.05 < K d . 

appear [33J. Therefore, it is very desirable to develop tools which quantify the degree 
of sticky motion, hyperbolicity and ergodicity. Following ideas from [6] the purpose of 
the present work is to quantify the degree of sticky motion in the phase space of higher- 
dimensional systems. To fulfill this task we have systematically analyzed small changes 
in the FTLE distributions using the skewness k^, which detects the asymmetry of the 
distribution, the kurtosis k 4 , which measures the flatness of the distribution and the 
normalized number of occurencies of the most probable finite time Lyapunov exponent 
P\, which indicates to which extent all trajectories started are under the influence of 
the same FTLE. 

Using coupled standard maps with nonlinearity parameter K we show that the 
quantities Ks(K), k^{K) and P\(K) are very sensitive to detect the sticky motion and 
capable to characterize the dynamics when the dimension d of the system increases 
(d = 2,4,6,10,20). Besides detecting accurately all cases of sticky motion, we found 
the following additional interesting results: i) (d = 6) the largest FTLE is most sensitive 
to the sticky motion; ii) (d = 10, 20) the smallest FTLE defines the critical nonlinearity 
parameter Kd, below which disconnected invariant regions in phase-space are expected 
and An\K) behaves nonlinearly. For K > K c a clear transition from the quasiregular 
to the totally chaotic dynamics occurs simultaneously in all unstable directions. In 
the interval Kd < K < K c all quantities detect the presence of sticky motion and 
distributions for different FTLEs behave similarly, meaning that the sticky motion 
affects all unstable directions proportionally by the same amount and a kind of common 
behavior of the distributions is observed. Physically what happens is that the linear 
momentum of all sites is constant at the sticky motion. 

Close-to-zero FTLEs play also an important role in other higher- dimensional 
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system, for example in connection to hydrodynamic modes [3U [35] (among others) and 
the analysis of hyperbolicity [HI ES] • Results of the present work contribute to categorize 
the rich and complicated quasiregular dynamics of higher-dimensional conservative 
systems and also to understand how nonlinearity is distributed along different unstable 
directions [T9| 137]. 
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